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We construct a model in which neutrinos and anti-neutrinos acquire the same mass but slightly 
different energy dispersion relations. Despite CPT violation, spin-statistics is preserved. We find 
that leptogenesis can be easily explained within this model, without upsetting the solar, atmospheric 
and reactor neutrino data. Leptogenesis occurs without lepton number violation and the non- 
equilibrium condition. We consider only three active Dirac neutrinos, and no new particles or 
symmetries are introduced. 
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Introduction. Despite the crucial significance of CPT 
symmetry in the conventional quantum field theory, it 
has been shown that string interactions may induce cou- 
phngs between Lorentz tensors and fermions in the low- 
energy 4D effective lagrangian [H. When the appropri- 
ate components of these Lorentz tensors acquire non-zero 
vacuum expectation values, they lead to a spontaneous 
CPT violation. 

Recently, the result from MINOS suggests a ten- 
sion between the oscillation parameters for and t'^ 
disappearance More precisely, at 90% confidence 



level, 
and 



it reports that: 
I Am^2 I = 3.36 



Am§2 



2.35 



-1-0.45 



+0.11 
0.08 



X 10-3 eV^ 

-0.40 ^ 10^3 eV"' , together with 
sin2( 26*23 ) > 0.91 and sin2(2^23) = 0.86 ± 0.11. This 
substantial difference in the neutrino and anti-neutrino 
mass-squared spUttings, if persists, may indicate CPT 
violation in the neutrino sector. 

The quest to explain MINOS could serve as a phe- 
nomcnological motivation for CPT violation in the neu- 
trino sector. In fact, there have been some works along 
this direction Q. But this motivation is not unique, as 
addition of sterile neutrinos Q or non-standard neutrino 
interactions Q may also provide an explanation. 

An earlier suggestion for CPT- violating neutrinos was 
due to the unresolved neutrino data in the Liquid Scintil- 
lator Neutrino Detector (LSND) @. In LSND, the flux 
of lies in the energy range 20MeV < E < 52.8 MeV, 
with an average about 40 MeV. The transition probabil- 



ity of 



has been measured. The data indicate a 



3.8 a excess of events, which, if interpreted as originat- 



ing from 



: oscillation, would imply a mass-squared 



splitting larger than 0.1 eV^. This is at odd with the data 
from SNO Q, KamLAND Q and Super-Kamiokande 0. 
It was then suggested that CPT-violating neutrinos may 
reconcile all the data [loj . 

In MiniBooNE, the e nerg y of neutrinos or anti- 
neutrino is of order GeV [11|. Recently, the — )• i^e 
search in MiniBooNE has found no evidence for an ex- 
cess of Ve in the — > i>e search [llj , and it reveals that 
a mass-squared spUtting smaller than 0.1 eV^ is allowed. 



However, their Vu 



study appears to be consistent 



with LSND [Ij] , and a mass-squared splitting of at least 
IQ-^eV^ is required. Again, this seems to suggest CPT 



violation in the neutrino sector. 

In this article, we would like to provide a new and 
simple model with only three CPT-violating active Dirac 
neutrinos, which explains leptogenesis without lepton- 
numbcr violation. At the same time, we show that this 
model is consistent with all of solar, atmospheric and 
reactor neutrino data. 

Neutrinos and CPT Violation. Parallel to the most 
popular framework advocated in [l3| . we propose a new 
type of Lorentz and CPT violations in the neutrino sec- 
tor. Our model is described by the following Lagrangian: 

L = Va{i5af)$- ruap ) V0 + 7° 3 A Ob vp , 

(1) 

where a,/? = e,fj,,T, map is the mass mixing matrix, 
\ap are dimensionless parameters characterizing kinetic 
mixings between different flavors of neutrinos, A^B = 
0,1,2,3 are spacetime indices and T^^ is a constant 
background tensor which breaks Lorentz invariance. For 
a given scalar function /, we deflne the unit operator Ba 
as the following: 9^/= (^o/, V/) = (^o/, ^) ■ 
Under Lorentz transformation, the unit operator trans- 
forms as Ba' ~ g^Ai Oa ■ The components of T^^ can 
be chosen such that T^^ = S"^^ ii A,B = 1,2,3 and 
= otherwise. In this way, ([1]) reduces to 



L 



map )vp+ Va Kp ^"V -Viyp. 



(2) 



As the neutrino held will be expanded as a linear combi- 
nation of the plane- wave solutions e*'^'^, the unit gra- 
dient operator V essentially operates as 



vv „± i p-x 

V e±^P-^ = = ±ipe^'P-^. 

I V e± 'P-^ I 



(3) 



Obviously, the operator V is ill-defined at zero- 
momentum p= 0. However, the composite operator V-V 
is well-defined because in the limit p — >■ 0, V • Ve^*^'^ 
vanishes. This is true regardless of the ordering of V 
and V, namely V • Ve^*^'^ also vanishes in the limit 
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We emphasize that the Lagrangian ([2]) is hermitian and 
renormahzable. The new operator with Xap 7° breaks C 
but preserves P and T, and so violates CPT. Since CPT 
violation implies Lorentz violation IJ] , this operator also 
breaks Lorentz invariance. While the particle Lorentz 
invariance is broken, the observer Lorentz invariance is 
preserved and so this new operator is consistent with the 
analysis provided by [3l ■ For a symmetric Xap , this op- 
erator will be identically zero if we consider Majorana 
neutrinos, because they do not have a vector current. 
If Xap is anti-symmetric, Majorana neutrinos will be al- 
lowed but CPT will no longer be violated. Since wc are 
interested in neutrino CPT violation, we arc essentially 
considering Dirac neutrinos in this article. 

We assume that the mass mixing and kinetic mixing 
matrices commute with each other, and so we can diago- 
nalize them simultaneously. Upon diagonalization by the 
usual unitary transformation, we obtain neutrino mass 
eigenstates and the Lagrangian becomes 



C = Va\i $ -ma) Sab + Va K Sab 7° V • V I/fc . 



(4) 



where a,b = 1,2,3. The corresponding energy disper- 
sion relations for the neutrino mass eigenstates are de- 
termined to be 



Ea = \/p'^ + + XaP 7 for neutrinos, (5) 
Ea = \/ p'^ + rUa — XaP , for anti-neutrinos. (6) 

We expect Aa ^ 1 to be consistent with current experi- 
ments. As a result, neutrino and anti-neutrinos acquire 
the same mass but slightly different energy dispersion 
relations. This is in contrast to the conventional sense 
of CPT violation in the neutrino sector, which requires 
neutrinos and anti-neutrinos to acquire different masses 

m- 

Since neutrino and anti-neutrinos acquire different en- 
ergy dispersion relations, the usual expansion of field op- 
erators in terms of creation and annihilation operators 
would have to be modified. The neutrino field operators 
are defined as 



(2^ 



E 



+ bl{p)v,{p) 



d^p 
(2^ 



E 



bsip)vsip) 



-tpx 



+ al{p) u,{p) 



(7) 



(8) 



where = Ep, = Ep and we have suppressed 
all the flavor or mass indices for generality. The cre- 
ation and annihilation operators can be imposed to obey 



the usual anticommutation relations: {ar{p), al{q)} = 
{&r(p), &|(q)} = (27r)3 5^,5(3)(p - q). This together 
with the usual sum rules for the spinors Us(p) and 
Us(p) lead to the equal-time anticommutation relations 
for the field operators: {i^(x), J^^(y)} = (^'■■^•'(x — y) and 
Mx),Ky)} = {^t(x), ^t(y)}^o. 

If we start from the Lagrangian ([2]) or Q and compute 
the conjugate momentum operator Tr(x) associated with 
the field operator, we obtain 7r(a;) = i i^^(a;). The canoni- 
cal quantization rule requires {i^(x), 7r(y)} =i(5'3^(x— y) 
and hence {i^(x), i'^(y)} — i5'-'^^ (x — y). This is obviously 
consistent with what wc derived from the neutrino field 
operators iy{x) and ^^{x) directly, and so the internal 
consistency of the entire construction is established. 
Therefore, we conclude that despite CPT violation in 
our model, spin-statistics is preserved. In fact, the 
above discussion reveals that any interaction term that 
breaks CPT but docs not contain dt v can preserve 
spin-statistics. 

Leptogenesis. A successful baryogenesis needs a pro- 
cess which satisfies all of the three Sakharov conditions 
[lit simultaneously: baryon number violation, C and CP 
violations, and non-equilibrium condition. One remark- 
able way to explain the observed baryon asymmetry in 
the universe is through leptogenesis. The main idea is 
that lepton asymmetry is preferentially generated in the 
very early universe. It is then partially transformed into 
baryon asymmetry by the sphaleron process [l7| which 
violates both baryon number (B) and lepton number (L). 
The analogous Sakharov conditions for leptogenesis are 
similar but with baryon number violation replaced by 
lepton number violation. In the standard paradigm of 
leptogenesis [18|, a heavy right-handed Majorana neu- 
trino decays into leptons and Higgs. This decay process 
is both L-violating and CP-violating. Interestingly, the 
right-handed Majorana neutrino is also responsible for 
explaining the smallness of neutrino masses through the 
see-saw mechanism figj . 

In contrast to the standard leptogenesis, CPT violation 
allows the Dirac left-handed neutrinos and right-handed 
anti-neutrinos to develop an asymmetry even at thermal 
equilibrium: 



dp 
2^ 



P 



1 



1 



(9) 



where we have set the Boltzmann constant fcs = 1 for 
convenience, T is the temperature, rii, and rip are the 
Fermi-Dirac distribution for neutrinos and anti-neutrinos 
respectively. Since spin-statistics is preserved in our 
model, we are safe to use the Fermi-Dirac distribution. 

With a given temperature T, the integrand in ([9]) is 
suppressed unless p ^ T. Thus, if ^/X^T ^ nia (which 
will be evidently justified in a moment), we can approxi- 
mate E Ki {1 + Xa)p and £" « (1 — Aa)p in the integrand. 
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C(3)t3 + 0(A3), (10) 



Performing the integration over p and keeping only the 
leading order, wc obtain the neutrino asymmetry 

Z TT^ 

for \/\iT ^ iria, with ^(3) ~ 1.202 being the Riemann 
zeta function. 

At the thermal equilibrium, the entropy per comoving 
volume is conserved. The entropy density is given by 
s (27rV45)5*r3 0. For T > 100 GcV, we have 

~ 106. Thus, the total neutrino asymmetry to entropy 
density ratio is 



E 



(11) 



A successful leptogenesis requires this ratio to be of order 
10^^'', which in turn requires 



10" 



(12) 



This is obviously valid if Ai = A2 = A3 = A which is 
implied from the lepton-number preserving case with a 
diagonal Aq,^ = A(5a^ in Thus, even if we keep a 
general and non-diagonal A^^ which violates the individ- 
ual neutrino lepton numbers, this violation is irrelevant 
for leptogenesis. We conclude that our model is capa- 
ble of generating the correct amount of lepton asymme- 
try without lepton number violation and non-equilibrium 
condition. 

As in the standard paradigm of leptogenesis, this pre- 
existing neutrino asymmetry will be partially converted 
into baryon asymmetry through the (B-l-L) -violating but 
(B-L)-preserving sphaleron processes (which are signifi- 
cant for T > 100 GeV). When the chemical equilibrium 
is reached, the baryon asymmetry is equal to [21| 



"^""^ = - 0.35 ^ 



a=l 



IQ- 



(13) 



As the sphaleron processes freeze out below the elec- 
trowcak scale, this baryon asymmetry will be perma- 
nently built into the quark sector. The quarks are con- 
fined to form baryons as the universe cools below the 
QCD phase transition scale (about 150 MeV), and this 
asymmetry becomes what we observe today. 

In retrospect, we confirm that since baryogenesis in our 
model occurs above the electroweak scale, the condition 
\/\^T ^ 771a and hence the approximations i? « (1 + 
Xa)p and E Ki {\ — \a) p are justified. 

We remark that if, on the contrary, one assumes 
left-handed neutrinos and right-handed anti-neutrinos to 
have different masses nia and ttTq, as was done by 
to resolve LSND, then the baryon asymmetry to entropy 
density ratio would have gone as 



10-^ E 



a=l 



rp2 



(14) 



For T > 100 GeV at which sphalerons are effective, we 
require — > (100 MeV)^ to ensure a successful 
baryogenesis, which is incompatible with the mass scale 
suggested by LSND or any other experiments. But we 
emphasize that our model of CPT-violating neutrinos 
predicts the correct amount of baryon asymmetry. 

In fact, an earlier idea of CPT-odd leptogenesis has 
been explored in 22 1 . The authors of (22j considered non- 
renormalizable dimension-5 operators (involving heavy 
Majorana neutrinos) that are CPT-violating and lepton- 
number violating. In comparison, our new idea of lepto- 
genesis from CPT violation is unique in the sense that 
both new particles and lepton number violation are not 
required. 

Furthermore, 



22 



listed a set of constraints on CPT- 
violating dimension-5 operators in the fermionic sector 
of the Standard Model H. In all of these works, CPT 



violation is achieved by the existence of a constant 
background vector. The constraints on the dimensionful 
coupling constants are thus derived. On the contrary, we 
are considering a renormalizable CPT-violating operator 
in the current paper. Now, CPT violation is achieved by 
the existence of a constant background tensor and the 
coupling constant is dimensionless. So it is not obvious 
that the constraints from [i^l are directly applicable 
to our work. We plan to find similar constraints in a 
forthcoming article. 

Implications for Neutrino Experiments. In the con- 
ventional neutrino oscillation formulae, the oscillation 
frequency is proportional to A.Eab = Ea — Eh, with 
a,b = 1,2,3. If Lorentz invariance and CPT are both 
preserved, the conventional energy dispersion holds, and 
the frequency oscillation is given by 



1 

2E 



ab I 



(15) 



where ^m^i, = rn^ — and E Ea ~ Eb because 
neutrinos are relativistic. However, in our model, neu- 
trinos and anti-neutrinos acquire the energy dispersions 
according to ((5)) and (jH) respectively. This means that 

AEab = ^[Aml, + 2E'{Xa-Xb)] , (16) 

AEab - ^[Aml,^2E'{Xa^Xb)] . (17) 

As a result, to confront our model with experiments, any 
experimental constraints on Am^f, will have to be re- 
interpreted as constraints on 

AM!,iE) ^ Ami, +2E'iXa-Xb), (18) 
AMl,iE) EE Ami, -2E\Xa-Xb). (19) 

Our model only modifies the usual energy dispersions 
of neutrinos and anti-neutrinos, but not the mixing an- 
gles. We adopt the usual mixing angles extracted from 
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solar (SNO), atmospheric (Super-Kamiokande) and re- 
actor (KamLAND, CHOOZ [2^) neutrino experiments. 
This means that we take sin^(20i2) ~ 0.8, sin^(26'23) ~ 
0.9 and sin2(26'i3) < 0.15. SNO and KamLAND have 
measured the survival probabilities of v^, — >■ Vf, and 
— > respectively, with both and v^, being in the 
MeV scale. Besides, Super-Kamiokande (SuperK) has 
measured the oscillation probability for atmospheric neu- 
trinos with energy from a few GeV up to 100 GeV (al- 
though the detector is not able to distinguish neutrinos 
from anti- neutrinos in the flux). 

By assuming the usual mixing angles, we are required 
to satisfy the following constraints from SNO and Kam- 
LAND respectively: 

AM|i ( MeV ) « 7.6 X 10"^ eV^ , (20) 
AM2i(McV) w 7.6 X 10^5 eV^ (21) 

Apparently, MINOS may be explained in the following 
way. In MINOS, the average neutrino or anti-neutrino 
energy is about GeV. So the constraints on the mass- 
squared splittings are translated into 

I AM|2( GeV ) I = 2.35 x 10"^ eV^ , (22) 
|AM32(GcV)| = 3.36 X 10^3 eV^ (23) 

For normal mass hierarchy (Am|2 > 0), the above con- 
ditions ([22]) and ^ can be satisfied if A?77|2 w 2.86 x 
10~'^eV^ and A2 — A3 « 2.5 x 10~^^. For inverted mass 
hierarchy (Ato§2 < 0)) these conditions are satisfied if 
Aml^ « -2.86 X IQ-^eV^ and A3 - A2 « 2.5 x 10-^2. 
The validity of using the 2-neutrino oscillation formula 
to perform the data fitting also requires AM|]^ ( GeV ) sa 

AM2i(GeV) « 7.6 x lO^^eV^ Thus, MINOS can be 
explained without upsetting the solar and reactor neu- 
trino data if 

Am^i « 7.6 X 10^5 eV^ ; | Am^2 I ~ 2.86 x 10"^ eV^ 
Ai = A2 ; |A3-A2|«2.5x IQ-^a . (24) 

However, the parameters in (|24p would imply 

I /^Ml^{ 100 GeV) I - I AM32( 100 GeV) | - eV^ (25) 

Since the SuperK data indicate a rather fiat oscillation 
spectrum up to the high energy region ('--^lOO GeV), the 
parameters in (|24p obviously predict too many oscilla- 
tions for high energy neutrinos or anti-neutrinos at Su- 
perK. Therefore, they are excluded by the SuperK data. 

The only way for our model to be consistent with all of 
the solar, atmospheric and reactor neutrino data is the 
simplest case with Ai = A2 = A3 = A. As mentioned 
earlier, this corresponds to the lepton-number preserving 
case with a diagonal Xap = XSa/s in ((2). In this case, 
leptogenesis is still explained, although we will have 

AM!,{E) = Ami = ^KbiE) = Ami, ■ (26) 



So the effect of CPT violation is completely invisible in 
all the neutrino oscillation experiments. The frequencies 
of neutrino oscillations predicted by our model and the 
conventional theory of neutrino oscillation are exactly 
the same. Of course, this would imply that our model 
cannot explain MINOS and the — De "anomaly" in 
LSND and MiniBooNE. 

Conclusion. We construct a new model in which 
neutrinos and anti-neutrinos acquire the same mass but 
slightly different energy dispersion relations. This sim- 
ple model of neutrino CPT violation explains leptogene- 
sis easily, without lepton number violation and the non- 
equilibrium condition. Also, it is consistent with all of 
the solar, atmospheric and reactor neutrino data. In ad- 
dition, according to FIGURE 13.10 in [25j, our model is 
also consistent with all other neutrino experiments such 
as KARMEN (40 MeV), Bugcy (MeV), CDHSW (GeV), 
NOMAD (50 GeV), Palo Verde (MeV), etc. 

It would be interesting to generalize the idea of the 
current model and see if MINOS and the i/^ — >■ De 
"anomaly" in both of LSND and MiniBooNE can be 
explained as well. We will explore this possibility in a 
forthcoming article. 
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